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Helium accreting CO white dwarfs witli rotation: helium novae 

instead of double detonation 
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Abstract. We present evolutionary models of helium accreting carbon-oxygen white dwarfs in which we include the effects of 
the spin-up of the accreting star induced by angular momentum accretion, rotationally induced chemical mixing and rotational 
energy dissipation. Initial masses of 0.6 M0 and 0.8 M0 and constant accretion rates of a few times 1O^*M0 /yr of helium rich 
matter have been considered, which is typical for the sub-Chandrasekhar mass progenitor scenario for Type la supernovae. It is 
found that the helium envelope in an accreting white dwarf is heated efficiently by friction in the differentially rotating spun-up 
layers. As a result, helium ignites much earlier and under much less degenerate conditions compared to the corresponding 
non-rotating case. Consequently, a helium detonation may be avoided, which questions the sub-Chandrasekhar mass progenitor 
scenario for Type la supernovae. We discuss implications of our results for the evolution of helium star plus white dwarf binary 
systems as possible progenitors of recurrent helium novae. 
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1. Introduction 

Type la supernovae are of key importance for the chemical evo- 
lution of galaxies, as they are a major producer of iron group 
elements (e.g. Nomoto et al. 1984; Renzini [19991 ). They were 
also found to be excellent distance indicator and have become 



an indispensable tool in cosmology (Phillips 119931 Hamuy et 
al. '1996' Branch ' 1998'!. The recent suggestion of a non-zero 
cosmological constant is partly based on observations of SNe 



la at high redshift (Leibundgut 1200 1 i . Given that distance de- 
terminations at high redshift through SNe la depend on the as- 
sumption of the homogeneity of SNe la light curves through- 
out the ages, an understanding of the possible diversity of their 
progenitors is crucial to evaluate this approach. Nevertheless, 
the progenitors of Type la supernovae have not been identi- 
fied yet, and the debate on their exact nature continues (e.g., 
LivioE^. 

One possibility to obtain a SN la is the detonation of the 
degenerate helium layer accumulated on top of a CO white 
dwarf due to mass transfer from its low mass helium star 
companion in a close binary system, which triggers a car- 
bon detonation in the white dwarf core. This is the so called 
double detonation or sub-Chandrasekhar mass scenario for 
SNe la, as it may allow to explode white dwarfs with masses 
well below the Chandrasekhar mass (e.g. Nomoto [T982I 
Fujimoto 1982 Limongi & Tornambe 1991, Livne 1990" 
Woosley & Weaver .1994. Livne & Arnett .1995J . While the 



capability of the helium detonation to ignite the CO core is 
still debated (e.g., Livio^OOl ), the helium detonation by itself 
would produce an explosion of supernova scale. 

Currently, the sub-Chandrasekhar mass scenario is not fa- 
vored as a major source of SNe la mainly because the light 
curves and spectra obtained from this model are not in good 



agreement with observations (e.g. Hoflich & Khokhlov ll9961 
Nugent et al. I1997t . Especially, the predicted presence of 
high velocity Ni and He is most stringently criticized (e.g. 



Hillebrandt & Niemever l2000l Livio l200lJ . 

On the other hand, stellar and binary evolution theory pre- 
dicts a realization frequency of binary systems such as helium 
star cataclysmics — which might produce double detonating 
sub-Chandrasekhar mass white dwarfs — which amounts to a 
few times 10"'^ yr"^ per galaxy (e.g. Iben & Tutukov ll99U 
Regos et al. 120021 which is comparable to the expected total 
SN la rate in the Milky Way. This raises the question why such 
explosions are practically never observed. 

We note that the sub-Chandrasekhar mass SN progenitor 
models which have been constructed so far neglected the ef- 
fects of rotation, which can be one of the primary factors 
determining the evolution of stars, in particular of massive 
stars (Langer 1998. Maeder & Mevnet 2000 l. Iben & Tutukov 
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J1991I I pointed out that rotation may indeed be important in he- 
lium star cataclysmic systems. Yoon & Langer f2002' '2004^ 
and Yoon et al. (2004i showed that effects of rotation might 
be essential for the evolution of accreting white dwarfs when 
the accreted matter contains a high specific angular momen- 
tum. The induced spin-up was found to change the white 
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dwarf structure significantly and to produce rotationally in- 
duced chemical mixing. 

In this paper, we suggest that rotation could play a key 
role in helium accreting white dwarfs such that in model which 
would produce a helium detonation this phenomenon is com- 
pletely avoided when the white dwarf spin-up is considered. 
After explaining the numerical method and physical assump- 
tions of the present study in Sect. |2l we investigate the evo- 
lution of helium accreting carbon-oxygen white dwarfs with 
accretion rates of ~ 10^^ Mq/jy, with the effects of rota- 
tion considered, in Sect. |3] Implications of our results for he- 
lium novae and neutron capture nucleosynthesis are discussed 
in Sect. 13 Our main conclusions are summarized in Sect.|5] 

2. Numerical method and physical assumptions 

We have computed the numerical models with a hydrodynamic 
stellar evolution code (Langer et al. ll988> . which incorporates 
the effect of the centrifugal force on the stellar structure and ro- 
tationally induced transport of angular momentum and chemi- 
cal species due to the dynamical and secular shear instability, 
the Goldreich-Schubert-Fricke instability and the Eddington- 
Sweet circulations (Heger et al. 2000, Yoon & Langer 2004 ). 
Conservation of angular momentum and energy of viscous flu- 
ids requires dissipation of rotational energy as angular momen- 
tum is transported by viscous friction in differentially rotating 
layers (e.g. Landau & Lifshitz [l984> . This effect is considered 
in our calculations following Mochkovitch & Livio ( I1989> as: 
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where uj is the angular velocity, r the radius and t'turb the turbu- 
lent viscosity due to the above mentioned rotationally induced 
instabilities (Heger et al. 120001 Yoon & Langer l2004b . 

Nuclear networks for the changes in chemical composi- 
tion and the nuclear energy generation include more than 60 
reactions (see Heger et al. 12000! for more details). In particu- 
lar, the i'*N(e-, z/)"C(a, 7)1^0 reaction (hereafter NCO reac- 
tion), which becomes active when p > 10^ g/cm^ (Hashimoto 
et al. I1984> . has been newly included for this study. We have 
used the ^^N(e~ , i^^^C reaction rate given by Martinez (2002, 
Private communication) and followed Caughlan and Fowler 



J1988> for the ^'*C(a, 7)^*0 reaction rate. The accretion in- 
duced heating is described following Neo et al. (1977), and the 
accreted matter is assumed to have the same entropy as that of 
the surface of the accreting star. 

Two initial masses, 0.6 and 0.8 M0, are considered for the 
CO white dwarf models. Since isolated white dwarfs are gen- 
erally found to rotate with a surface velocity of Vg ^ 40 km/s 
(Heber et al.T^, Koester et al.lW Kawaler'^UU^^, the ini- 
tial rotation velocity of our models is assumed to be as slow 
as 10 km/s at the white dwarf equator (see also Langer et 
al. '19991. Other physical properties of the white dwarf initial 
models are summarized in Tabled While most of our simula- 
tions start with a cold white dwarf with log Ls/L© ~ —2.0, an 
initially hot white dwarf with logLg/L© ~ 2.508 is also con- 
sidered for one model sequence (TC, Tabled cf. Sect. 14. 1> . The 



Table 1. Physical quantities of the initial white dwarf models: 
mass, surface luminosity, central temperature, central density, 
radius and rotation velocity. The hot white dwarf model in the 
third row is only used for sequence TC (see Table|2}. 



A/wD.init 


log is,init/LQ 


Tc,init 


Pc,init 


-RwD,init 


^rot ,init 


Mq 




10^ K 


10*^ g/cm^ 


R© 


km/s 


0.6 


-2.049 


1.69 


3.64 


0.0126 


10 


0.8 


-2.024 


1.59 


10.7 


0.0101 


10 


0.8 


2.508 


14.4 


6.44 


0.0186 


10 



accreted matter, received with two different constant accretion 
rates of M = 2 x 10"^ and 3 x 10"^ Mo/yr, is assumed to 
have Y = 0.982 and Xn = 0.012, where Y and Xn are the 
mass fraction of helium and nitrogen, respectively. 

In a close binary system, the white dwarf is believed to re- 
ceive matter through a Keplerian accretion disk if its magnetic 
field is not considerable. The accreted matter may thus carry 
an amount of specific angular momentum which corresponds 
to critical rotation at the white dwarf equator. However, con- 
tinuous angular momentum gain under these conditions leads 
to over-critical rotation soon after the onset of mass transfer 
in the outer part of the accreting star (Yoon & Langer 2002}. 
Therefore, we limit the angular momentum gain such that the 
accreting star may not rotate over-critically, as follows: 



Ja 



/ • jKoplcr if Ws < / • WKeple 
if Ws = / • WKeple 



(2) 



where denotes the surface velocity at the white dwarf equa- 
tor, jacc the specific angular momentum of the accreted matter, 
and WKcpicr and jKcp the Keplerian value of the surface velocity 
and the specific angular momentum at the white dwarf equa- 
tor, respectively. The dimensionless parameter / represents the 
fraction of the Keplerian value of the specific angular momen- 
tum which is contained in the accreted matter, and also the 
maximum fraction of the critical rotational velocity with which 
the white dwarf equator may rotate. 

The use of / = 1 might be the most natural choice to de- 



scribe the realistic situation, as Paczynski (11991> and Popham 
& Narayan (1991) argue that a critically rotating star may con- 
tinue to accrete matter from a Keplerian disk by transporting 
angular momentum from the star back into the disk due to tur- 
bulent viscosity. However, as the correct treatment of close- 
to-critical rotation is beyond the capabilities of our numerical 
code, we also consider the case / < 1. As discussed in Yoon & 
Langer (2004), the correction factors of /p and /t which are 
included in the stellar structure equations for describing the ef- 
fects of rotation (cf. Heger et al. 12000> are limited to 0.75 and 
0.95, respectively. This limit corresponds to a rotation rate of 
about 60% of critical rotation, up to which our one dimensional 
approximation in computing the effective gravitation potential 
can accurately describe the structure of the rotating star In our 
models with / = 1, where the outer envelope rotates close to 
critically, the centrifugal force is accordingly underestimated. 
Although the region which rotates faster than 60% critical con- 
tains only httle mass (see Sect. |3j, we also consider the case 
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Table 2. Properties of the computed model sequences. MwD.init ■ initial mass, M : accretion rate, / : fraction of the Keplerian 
value of the accreted specific angular momentum (see Eq. |2}. £^diss ■ Rotational energy dissipation due to frictional heating. 
Yes (or No) means that this effect is considered (or not). AMhc : accumulated helium mass until the helium ignition point. Tc 
and pc : central temperature and density in the last computed model. Thc : maximum temperature in the helium envelope in 
the last computed model. pHc and 77H0 : density and degeneracy parameter at the position of the maximum temperature in the 
helium envelope in the last computed model. The last column indicates whether the model sequence will finally result in helium 
detonation or not. 



No. 


-MwD.init 


M 


/ 


-Bdiss 


AMhs 


Tc 


Pc 




pHe 


7?Ho 


Detonation? 




Mq 


iU MQ/yr 






M0 


iU Jv 


iU g/cm 


1 nS 
iU K 


iU g/cm 






Nl 


0.6 


2.0 






0.229 


0.60 


10.9 


2.0 


1.48 


11 


Yes 


IN/ 


U.D 








U. ioj 


U.D4 


0.0 


1 A 

1.4 


V.JJ 


y 


No 


IN J 


C\ S 
u.o 


z,u 






U. iOo 


U.J J 


ZD. 1 


z.u 


1 

1 .Jo 




lea 


N4 


0.8 


3.0 






0.150 


0.46 


23.6 


1.2 


1.50 


19 


Yes 


Rl 


0.6 


2.0 


1.0 


Yes 


0.080 


0.39 


4.4 


1.2 


0.22 


6 


No 


R2 


0.6 


3.0 


1.0 


Yes 


0.061 


0.23 


4.2 


1.4 


0.15 


4 


No 


R3 


0.8 


2.0 


1.0 


Yes 


0.022 


0.14 


11.0 


1.5 


0.15 


4 


No 


R4 


0.8 


3.0 


1.0 


Yes 


0.017 


0.14 


10.6 


3.8 


0.03 


0.3 


No 


TAl 


0.6 


2.0 


0.6 


Yes 


0.149 


0.50 


5.97 


1.2 


0.42 


9 


No 


TA2 


0.6 


3.0 


0.6 


Yes 


0.106 


0.45 


4.70 


3.5 


0.16 


2 


No 


TA3 


0.8 


2.0 


0.6 


Yes 


0.073 


0.21 


14.2 


1.3 


0.38 


8 


No 


TA4 


0.8 


3.0 


0.6 


Yes 


0.045 


0.18 


12.7 


1.4 


0.28 


6 


No 


TBI 


0.6 


2.0 


1.0 


No 


0.402 


0.59 


9.6 


1.1 


1.14 


22 


Yes 


TB3 


0.8 


2.0 


1.0 


No 


0.233 


0.65 


21.3 


1.2 


1.36 


21 


Yes 


TC 


0.8 


1.0 


1.0 


Yes 


0.002 


0.95 


8.5 


2.8 


0.001 


0.06 


No 



/ = 0.6, apart from the case / = 1.0. With / = 0.6, the white 
dwarf never rotates faster than 60 % of the Keplerian value, and 
thus the stellar structure is accurately described throughout the 
white dwarf interior 

For comparison, we also compute rotating models where 
the rotational energy dissipation is neglected, as well as non- 
rotating models, with otherwise identical initial conditions. In 
Table 12 we list all computed model sequences. The index N 
in the model number indicates a non-rotating model sequence, 
while R denotes a rotating one with / — 1.0. The index TA 
is for rotating test models with / = 0.6, while TB is used for 
rotating test models with / = 1.0 without rotational energy 
dissipation (i.e., ediss.r = 0.0). The sequence TC designates 
the only model starting with a hot white dwarf. All sequences 
are computed up to the point where the accumulated helium 
shell ignites. 

3. Results 

We summarize the results of our simulations in Tabled Here, 
Thc, PHo and t^hc denote, for the last computed model, the 
maximum temperature in the helium envelope and the corre- 
sponding density and the degeneracy parameter (ip/kT, e.g. 
Kippenhahn & Weigert 1990). AA/hc gives the accumulated 
helium mass until helium ignition. 

At helium ignition we stop our calculations, and our models 
are per se not able to predict whether the helium burning devel- 
ops into a helium detonation or not. However, from the litera- 
ture (e.g., Woosley & Weaver [T994> we conclude that a helium 
detonation can not develop if the ignition density is smaller 
than ~ 10^ g/cm'^, due to the quenching of the thermonuclear 
runaway by expansion. For pHe < 10® g/cm'^, therefore, the 



helium ignition may result only in a nova-like shell flash, which 
will not be able to trigger core carbon ignition. 

Our results for the non-rotating model sequences are 
found in good agreement with similar models computed by 
Woosley & Weaver (1994) and Piersanti et al. (2001). For 
instance, for a sequence with the same initial mass and ac- 
cretion rate as for our sequence Nl, Piersanti et al. obtained 
AMhc = 0.244 Mq, which does not differ much from our re- 
sult of AMhc = 0.229 Mq. The slightly smaller value of the 
present study may be attributed to a small difference in the ini- 
tial nitrogen mass fraction, which triggers the NCO reaction 
when Pc ^ 10® g/cm^ and initiates helium burning. 

Fig. [2 illustrates the evolution of our white dwarf models 
with MwD.init = 0.8 Mq and M = 2 X 10"^ MQ/yr (sequences 
N3, R3 and TB3). In the non-rotating case, the temperature 
of the helium envelope continues to increase due to accretion 
induced heating. When the white dwarf mass reaches about 
0.9 Mq, the density at the bottom of the helium envelope starts 
exceeding 10® g/cm^ and the NCO reaction becomes active, 
accelerating the temperature increase. Finally, helium burning 
starts when A/wd reaches 0.968 Mq. The density at the he- 
lium ignition point is about 1.6 x 10® g/cm^ and a detonation 
is likely to follow. 

The evolution of the corresponding rotating model (R3) in 
Fig lDlooks similar, but helium ignites much earlier than in the 
non-rotating case. The maximum temperature in the helium en- 
velope reaches 10^ K when Mwd — 0.818 Mq. Helium burn- 
ing develops quickly thereafter, and the nuclear energy genera- 
tion amounts to logLHo/-^© — 5.0 when AfwD — 0.822 Mq. 
Only 0.022 Mq of helium has been accumulated by then, 
which is about 10 times less than in the non-rotating case. 
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Fig. 1. Evolution of accreting white dwarf models of sequence 
N3, R3 & TB3, for which MwD,init= 0.8 M© and A/ = 2 x 
10~® Mo/yr, in the log/?— log T plane. The boundary between 
the CO core and the helium envelope, defined though Y = 0.1, 
is indicated by a filled circle for each model. 

Importantly, the ignition density pHc is also found about 10 
times lower (1.5 x 10^ g/cm^) than in the non-rotating case, 
implying that a helium detonation may be avoided. 

The reason for the early helium ignition in the rotating se- 
quence is as follows. As explained above, the white dwarf ob- 
tains angular momentum carried by the accreted matter. The 



h . . . / I . ... . I ^. . / . . I J , , H 
0.65 0.70 0.75 0.80 0.85 

Mr/Mo 

Fig. 2. Upper panel: Angular velocity in units of the local 
Keplerian value as function of the mass coordinate in model 
sequence R3 when Mwd = 0.805, 0.812 and 0.821 Mq. The 
filled circles indicate the boundary of the CO core and helium 
envelope, as defined in Fig. [2 where the helium mass fraction 
is 0. 1 . Lower panel : Rotational energy dissipation due to fric- 
tion (Eq. Q at the same evolutionary epochs as in the upper 
panel. The thin solid line denotes the energy generation rate 
due to nuclear reactions when A/wd =0.821 Mq. 



angular momentum is thus transported from the surface into 
the interior by various rotationally induced hydrodynamic in- 
stabilities. In our white dwarf models, Eddington Sweet circu- 
lation, the secular shear and the GSF instability dominate in 
the non-degenerate envelope, while the dynamical shear insta- 
bihty is most important in the degenerate interior, as discussed 
in detail by Yoon & Langer f2004'i. The white dwarf interior 
is spun up progressively with time, giving rise to rotational en- 
ergy dissipation as shown in Fig. |2l In the upper panel of the 
figure, the angular velocity reveals a strong degree of differen- 
tial rotation. The consequent energy dissipation rate is given in 
the lower panel in the same figure. The total energy dissipa- 
tion rate integrated over the spun-up layers in the given models 
amounts to 8.8 Lq, 6.2 Lq, and 3.6 L© for AfwD = 0.805 Mq, 
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0.812 Mq, and 0.821 M0 respectively. Together with the ac- 
cretion induced heating, this additional energy supply speeds 
up the temperature increase in the helium envelope, leading to 
the earlier helium ignition compared to the non-rotating case. 

For comparison, the white dwarf evolution in sequence 
TB3, where the rotational energy dissipation is neglected while 
everything else is as in sequence R3, is shown in the third panel 
of Fig.n In this case, helium ignition occurs even later than in 
the corresponding non-rotating case (N3). This is due to the lift- 
ing effect of the centrifugal force, which reduces the accretion 
induced heating. The density at the bottom of the helium enve- 
lope when the helium burning is induced by the NCO reaction 
is about 1.36 x 10^ g/cm^, implying that a helium detonation 
might be triggered in this case. 

As shown in Fig. |2] the surface layers of the white dwarf 
models in sequence R3, where / = 1 is used, rotate faster 
than 60 % critical, above which the effects of rotation are un- 
derestimated by our computational method. In the models with 
MwD,init = 0.8 M0, i.e., in sequences R3 and R4, however, this 
fast rotating region is limited to only 1% of the white dwarf 
mass, which is not likely to affect the result significantly. On 
the other hand, about 7% in mass, which include a significant 
fraction of the helium envelope, exceed the limit at helium ig- 
nition in model sequences Rl and R2, for which MwD,init = 
0.6 M0 and f — I are adopted. Nevertheless, results from the 
corresponding sequences with / = 0.6 (TAl, TA2, TA3 and 
TA4), where the white dwarf models are forced to rotate below 
60% critical throughout the white dwarf interior, lead us to the 
same conclusion as in case of / = 1: helium ignites at such a 
low density (even though it is somewhat higher than in case of 
/ = 1) that a supernova event is unlikely to occur For instance, 
in model sequence TAl, we have pHe = 4.2 x 10^ g/cm^ 
when the maximum temperature reaches 1.2 x 10^ K. This is 
about 3.5 times lower than in the corresponding non-rotating 
case (phc = 1-48 x 10^ g/cm'^) and about two times higher 
than in case of / = 1 (pHo = 2.2 x 10^ g/cm^). 

4. Discussion 

4. 1. Connection to recurrent helium novae 

As discussed by Iben & Tutukov (1991 1 and Limongi & 
Tornambe (1991 ) , a relatively low helium accretion rate of 
~ 10^^ Mq/yt onto a CO white dwarf can be reaUzed, for 
example, in binary systems which consist of a 0.6 • • • 1.0 Mq 
CO white dwarf and a less massive helium star Mass trans- 
fer in such binary systems is driven by gravitational wave ra- 
diation, and the resulting mass transfer rate is found to be 
insensitive to the exact masses of white dwarf and helium 
star, at a few times 10"^ Mq/jt (Iben & Tutukov [T99T1 
Limongi & Tornambe 19911. Many of such binary systems, 
if not all, could produce a supernova which will appear as 
Type la (Taam [Mol Nom oto IT982I Iben & Tutuko v IT99T1 
Limongi & Tornambe' 1991' Woosley & Weaver' 1994' Livne & 
Arnett 1995), while the resulting light curves and spectra may 
be abnormal (Hoflich & Tutukov IT9961 Nugent et al. [T997t . 
As mentioned in the introduction of this paper, binary popu- 
lation synthesis models show that the production rate of such 



events may be comparable to the observed SN la rate (e.g. Iben 
& Tutukov [l9911 Regos et al. 2002). and therefore the appar- 
ent absence of such a peculiar type of supernovae has been a 
puzzling matter 

The results of the present study offer a possible solution to 
this problem: helium accreting CO white dwarfs with an ac- 
cretion rate of ^ 1G~® Mq/jt may not result in a supernova 
at all, but may instead produce nova-like explosions. The ac- 
cumulated helium mass of AMhc — 0.02 ■ • • 0.1 Mq in our 
rotating models indicates that in a low mass helium star + CO 
white dwarf binary system, nova explosions will occur recur- 
rently, with a period of about 10^ yr. 

As pointed out by Iben & Tutukov ( I1991> . various factors 
need to be considered for predicting the further evolution of 
such a binary system after the first nova outburst. First, any 
mass loss induced by the helium flash will affect the binary or- 
bit. I.e., the nova induced mass loss will widen the orbit and in- 
terrupt the mass transfer for some time (Iben & Tutukov ri991> . 
However, the angular momentum loss due to gravitational wave 
radiation will lead the helium star to fill its Roche lobe soon 
again (after 10^ yr) as discussed below. 

Second, the first nova outburst may heat up the white dwarf 
significantly, and the white dwarf may still be hot when the 
second mass transfer starts, compared to the case of the first 
mass transfer Model sequence TC, where an initially hot white 
dwarf with logLs/L© = 2.508 is adopted, was computed in 
order to investigate the effect of a pre-heated white dwarf. A 
value of M = 10~* Mq/yt for this sequence, which is smaller 
than in other model sequences, was chosen to consider the de- 
crease of the mass transfer rate due to the change of the or- 
bit as the binary system loses mass via helium nova outbursts 
(cf. Fig |3}. A comparison of AA/ho in sequence TC (AMhb 
~ 0.002 M0) with the accumulated helium masses in the se- 
quences R3 and R4 (AMue ^ 0.02 Mq) implies that AMho 
decreases by more than a factor of 10 if the white dwarf is pre- 
heated. Therefore, the second and any further helium flash may 
be much weaker than the first one. 

In an attempt to make more detailed conjectures on the evo- 
lution and final fate of binary systems of the considered kind, 
we made a simple experiment as follows. We constructed a bi- 
nary star model consisting of a zero-age helium main sequence 
star of 0.6 Mq and a CO white dwarf of 0.8 Mq in an 1.08 h 
orbit (Aorbit,init — 0.6 Rq, cf. Limongl & Tornambe II 99 U . 
Here, the white dwarf is approximated by a point mass. Since it 
appears likely that the accreted helium will be ejected from the 
system by the violent helium shell flash when about 0.02 Mq 
of helium is accumulated as implied by the results of sequences 
R3 and R4, we assume this for the evolution of our binary sys- 
tem. As discussed above, the subsequent helium flashes may 
occur with a smaller AA/ho- Therefore, for the subsequent evo- 
lution the white dwarf is assumed to lose the accumulated mat- 
ter at every time when AA/ho = 0.002 Mq is achieved. For 
comparison, a second binary evolution model with the same 
initial condition is calculated, with the assumption of no mass 
loss due to helium flashes. The evolution of the helium star and 
the change of the binary orbit due to the mass transfer, stellar 
wind mass loss and gravitational wave radiation are followed 
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Fig. 3. Mass transfer rate as function of the helium star mass, 
for a binary system consisting initially of an 0.6 Mq helium 
star and an 0.8 M© CO white dwarf (the latter here being ap- 
proximated by a point mass). Mass transfer is initiated by gravi- 
tational wave radiation. The solid line illustrates the case where 
mass loss due to helium flashes is considered, while the dashed 
line was computed assuming that no mass loss occurs. 
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Fig. 4. Evolution in the HR diagram of the helium star (initial 
mass 0.6 Mq) component of the binary system considered in 
in Fig.|3j for case with nova induced mass loss. 

by using a binary stellar evolution code (see Langer et al. 120001 
for more details about the code). 

Fig- El shows the evolution of the mass transfer rate in the 
considered system, as function of the helium star mass (A/hcs)- 
The mass transfer from the helium star starts when the helium 
mass fraction in the helium star center equals 0.39. The orbital 
period at this moment is about 39 min. The mass transfer rate 



initially increases to 3 x 10 * Mq/ji and has decreased to 
2.3 X 10"^ Mg/yr by the time AMho reaches 0.02 Mq. As 
the white dwarf is assumed to lose 0.02 Mq at this point due to 
the helium nova flash, the system becomes detached for about 
5.5 X 10^ yr, after which the helium star again fills its Roche 
lobe. Further-on, the white dwarf loses mass whenever AMne 
reaches 0.002 Mq, and thus the mass transfer is switched on 
and off repeatedly every 1.5 x 10'^ yr. We follow the evo- 
lution of the system until the helium star mass decreases to 

0.30 Mq. 

Fig. El shows the evolution of the helium star in the HR 
diagram. The helium star luminosity decreases continuously as 
it loses mass, and it will finally evolve into a white dwarf. At 
this stage, the binary system will resemble an AM CVn system. 

At the end of the calculation (i.e., A/hcS = 0.30 Mq), the 
central helium abundance has decreased to 0. 16, and more than 
0. 1 Mq of helium is still available in the envelope for further 
mass transfer. Since about 135 nova outbursts occurred until 
the end of the calculation, more than 180 recurrent nova out- 
bursts in total are expected to occur throughout the evolution 
of the considered binary system. This implies that heUum nova 
explosions in low mass helium star + CO white dwarf binary 
systems could be realized with a frequency of ^ 0.1 yr~^ 
in our Galaxy, given that such binary systems are being pro- 
duced at a frequency of ~ 0.001 yr~^ (Iben & Tutukov [l991l 
Regos et al. 2002 ). The recently discovered outburst of V445 
Puppis, which has been attributed to a helium nova (Ashok & 
Banerjee '2003' Kato & Hachisu '2003"l, may be a promising 
observational counterpart of such an event. The high observed 
carbon abundance in this system (Ashok & Banerjee I2003> 
might be explained by the rotationally induced chemical mix- 
ing in the accreting white dwarf, the effects of which are dis- 
cussed in the next section. 

4.2. Rotationally induced chemical mixing and 
neutron capture nucleosynthesis 

As already suggested by Iben & Tutukov (I1991> . strong helium 
flashes in accreting CO white dwarfs may activate neutron cap- 
ture nucleosynthesis, since neutrons can be provided via the 
^^Ne(a, 7i)^^Mg reaction. We note that favorable conditions 
for neutron capture nucleosynthesis might be achieved with ro- 
tation compared to the non-rotating case, due to the rotationally 
induced chemical mixing (cf. Langer et al. [T999i . 

Fig-Et gives the mass fraction of "'He, ^^C and ^^Ne in the 
initial model (MyvD = 0.8 Mq) of sequence R3, as well as in 
the model with AfwD = 0.822 Mq, which is immediate before 
the helium ignition. It is found that '^^C and ^^Nc have been 
significantly dredged up into the helium envelope when helium 
ignites, due to rotationally induced mixing. For instance, com- 
paring the two models shown in Fig. |5t reveals that the mass 
fraction of ^^C and ^^Ne at ~ 0.8 Mq has increased by 
about a factor of 10 when helium ignites. 

The ignition of helium induces convection. Fig. EJ) shows 
the chemical structure when Thc reaches 1.5 x 10^ K. At this 
point, the convective layer has a mass of AMconv — 0.01 Mq, 
extending from = 0.796 Mq to = 0.806 Mq. The 
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E (a) R3 (Minit = 0.8 Mo) 




0.76 0.77 0.78 0.79 0.80 0.81 0.82 

Mr/Me 

Fig. 5. (a) Chemical profiles in the outer layers of the white 
dwarf models of sequence R3 when Mwd = 0.8 Mq (initial 
model, thin lines), and AfyvD = 0.822 Mq (thick lines) which 
is immediate before the helium ignition. The solid and dotted- 
dashed lines give the mass fraction of ^He and ^^C respec- 
tively, as function of the mass coordinate. The dashed line de- 
notes the mass fraction of ^^Ne. (b) Same as in (a) but after the 
ignition of helium. The dotted line denotes the neutron density, 
for which the scale is given at the right side. 



abundance of the dredged up core material in the convective 
layer amounts to 0.15 for ^^C and 5.5 x 10^"^ for ^^Ne, respec- 
tively. Neutrons are released by the ^^Ne(a, n)^^Mg reaction, 
and a maximum neutron density of ~ 10^ cm"'^ is achieved at 
Mr = 0.76 Mq in our last model. This neutron density may 
be somewhat overestimated in our calculation, since we did 
not include the reaction ^'*N(n, p) ^"'C in our network (cf Siess, 
Goriely & Langer 2003i) . However, we note that alpha parti- 
cles are mixed significantly into the ^^Ne-rich region, and that 
large amount of neutrons will be released in the further evolu- 
tion of the flash which is not covered any more by our model, 
as hydrodynamic effects in the exploding layers may become 
important. Therefore, we leave the exploration of the neutron 
capture nucleosynthesis effects in our models for the future. 



5. Concluding remarks 

We have shown that the effects of rotation in helium accreting 
white dwarfs may be incompatible with the scenario of double 
detonations in sub-Chandrasekhar mass CO white dwarfs as 
possible SNe la progenitors. In helium accreting white dwarfs, 
we find the thermal evolution to be affected by viscous heat- 
ing due to differential rotation in the spun-up layers, such that 
helium ignition is induced at too low densities to develop a det- 
onation. This may give a plausible solution to the long stand- 
ing problem of the missing observational counterparts of sub- 
Chandrasekhar explosions, which are predicted to occur with a 
frequency comparable to the observed SN la rate. 

We discussed that binary systems consisting of a CO white 
dwarf and a less massive helium star may be possible progen- 
itors of recurrent helium novae (Iben & Tutukov 1991), which 
may be analogous to V445 Puppis (Ashok & Baneriee 120031 
Kato & Hachisu r2003t . After the first strong helium nova flash 
in such binary systems, rather mild nova outbursts are expected 
to occur recurrently with a period of ~ 10^ yr. The realization 
frequency of such a helium nova may be as high as ^ 0.1 yr~^ 
in our Galaxy. 

Rotation induces chemical mixing of ^^Nc and *He at the 
bottom of the helium envelope, which may provide interesting 
conditions for neutron capture nucleosynthesis to occur during 
the helium nova flashes. 

Finally, we note that other important mechanisms for the 
angular momentum redistribution in white dwarfs may exist 
than those considered in the present study. In particular, we ne- 
glected the possible role of magnetic fields, which may increase 
the efficiency of the angular momentum transport significantly 
(cf. Heger et al. 120031 Maeder & Meynet 2003 1. If the spin- 
up time scale is shorter than considered in the present study, 
the resulting shear strength will be weaker and the effect of ro- 
tational energy dissipation may not be as important as shown 
here. However, studies of magnetic effects in accreting white 
dwarfs have to be left for future investigations. 
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